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Abstract 

The origin of supermassive black holes in the galactic nuclei is quite 
uncertain in spite of extensive set of observational data. We review the 
known scenarios of galactic and cosmological formation of supermassive 
black holes. The common drawback of galactic scenarios is a lack of 
time and shortage of matter supply for building the supermassive black 
holes in all galaxies by means of accretion and merging. The cosmological 
scenarios are only fragmentarily developed but propose and pretend to an 
universal formation mechanism for supermassive black holes. 



Contents 

1 Introduction 

2 Theoretical fundamentals 

3 Observational signatures of SBHs 

3.1 Active galactic nuclei and quasars .... 

3.2 Dynamics of gas and stars near SBH . . . 

3.3 Galaxy — SBH correlations 

3.4 Black holes of intermediate mass 

3.5 Prospects for gravitational waves detection 

3.6 Astrophysical constraints on SBHs .... 

3.7 Dark matter spike in the Galactic center . 

* dokuchaev@inr.npd.ac.ru 
^ erosh@inr.npd.ac.ru 
■f sergeirubin@list . ru 



4 Astrophysical scenarios of SBH formation 

4.1 Dynamical evolution of galactic nuclei 

4.1.1 Non-dissipative contraction of star cluster 

4.1.2 Dissipative contraction of star cluster . . 

4.2 Collapse of supermassive star 

4.3 Collapse of neutron star cluster 

4.4 Merging of stellar-mass black holes 

4.5 Merging trees of galaxies and SBHs 



Cosmological scenarios of SBH formation 

5.1 Population III stars or supermassive stars 

5.2 Collapse of massive primordial clouds 

5.3 Pre-galactic SBH 

5.4 Adiabatic fluctuations at radiation-domination stage 

5.5 Clusters of BHs from closed domain walls 

5.5.1 General idea 

5.5.2 SBH formation in the Hybrid Inflation 

5.5.3 Evolution of massive cluster of primordial BHs |28 



6 Conclusions and Discussions 



1 Introduction 

The physics of Black Holes (BHs) is the most developed branch of general rela- 
tivity. There were elaborated numerous mechanisms for the formation of BH of 
different mass-scales and proposed a list of convincing observational signatures 
in favor of BH existence. A possible range of BH masses in the Universe is 
thought to be extremely wide: from the Plank mass, Mp\ ^ 10^^ g, up to the 
huge mass of the most luminous quasars, A/bh 10^°Mq. Conventionally we 
will distinguish the stellar-mass BHs with mass in the range ~ 3 — 100 Mq, the 
Intermediate Mass Black Holes (IMBHs) with mass in the range ^ 10"^ — IO^Mq 
and the Supermassive Black Holes (SBHs) with mass > IO^Mq. The stellar- 
mass BHs are rather definitely originated from the collapse of exhausted massive 
stars during supernova explosions. At the same time the origin of other types 
of BHs nowadays is quite uncertain. 

Modern astrophysics provides strong indications of the existence of BHs as 
among the remnants of evolved stars and in the centers of spiral and elliptic 
galaxies. Meanwhile to declare the real discovery of BHs we need the proof that 
the suspicious candidates have a distinctive BH feature — an event horizon. 

In contrast to the evident origin of stellar mass BHs the origin of SBHs with 
a mass exceeding IO^Mq in the galactic nuclei is quite unclear. With a present 
state of art of observations and theory we have nowadays only the list of more 
or less elaborated scenarios leading to the SBH formation. We will review here 
these known ideas and models with discussions of possible crucial signatures to 
discriminate between them and find out the universal mechanism (or several 
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mechanisms ) of SBH formation. The recent advance and current observational 
status of SBH searches is extensively reviewed e.g. in [H [5] . 

The SBHs is not those things that can be "touched by hands" or be in- 
vestigated thoroughly in the laboratory (it must be noted that formation of 
short-lived microscopic BHs is forecasted in particles collisions at LHC in ac- 
cordance with some extra-dimensional theories). Nevertheless, a supposition of 
their existence permits to explain a lot of observations. The Seyfert and radio 
galaxies, blazars, quasars and huge jets emitted from the galactic centers could 
be explained as various aspects of the same unique phenomenon: an activity of 
SBH in the galactic center. It is an interaction of SBHs and matter in the host 
galaxies that produces and explains different forms of Active Galactic Nuclei 
(AGN) phenomenon [Tl[l[3]. 

A recent discovery of very distant quasars reveals the deficiencies in the 
understanding of SBH formation. The SBHs are inevitably the main engine of 
quasars and the same fact of observation of so distant SBHs is a source of many 
questions. Indeed, how could SBHs as massive as 10^ — IQ-'^^M© be formed 
so early, at 2 = 6 — 7, i. e. at 800 million years after the Big Bang? After 
nucleation, a BH duplicates its mass during the time [6] 

tacc ^ 4 10^ yrs, 

if an accretion near the Eddington limit takes place. In this case the smallest 
period of time necessary to increase an initial BH mass Mq up to some value M 
is 

4 • 10^ , M 
t ~ In yrs. 

In 2 Mo 

Hence about billion years is needed to increase the BH mass from e. g. the 
Solar mass to the observable values. There is no room for all stages preceding 
nucleation of the Solar mass BHs — from a star formation to its explosion. 

Nevertheless, the multi-scale simulations indicate that luminous quasars at 
z ~ 6.5 could be formed within the framework of standard ACDM paradigm. It 
was shown [7] that the most massive halos after several mergers could reproduce 
in general the observed properties of SDSS J1148+5251, the most distant quasar 
detected at z ^ 6.42, [5]. 



2 Theoretical fundamentals 

The physics of BH is based on the Einstcin-Hilbert equations of general relativ- 
ity. Specific details and features of BH are investigated since 1916, when Karl 
Schwarzschild found the famous static spherically symmetric solution for a gen- 
eral relativistic body 9_. In spite of complexity of the subject and permanent 
disputes on the BH existence, scientific community has came to agreement on 
the general BH features. There are a few important formulas of BH physics 
which are useful in interpretation of observational data and astrophysical appli- 
cations. 
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A size of BH is characterized by its gravitational radius 

2GM ^ M 

., = ^.310^— err., (1) 

where G is the gravitational constant and M is a BH mass. Another important 
scale is the minimal radius Tgt of stationary circular orbit of a massive body 
moving around the BH: 

Tst = 3rg (2) 

in the case of nonrotating BH. 

The SBHs are the most bright objects in the Universe due to radiation of 
accreting matter. In a dense environment of galactic nucleus the mass of SBH 
is increased in several orders of magnitude during the lifetime of the galaxy. In 
the stationary spherically symmetric case (Bondi accretion) the growth of BH 
is defined mostly by the density of surrounded gas. The corresponding rate for 
BH mass growth is 

dMB G^APp 



dt 



Ana (3) 



where a dimensionless parameter a depends on the equation of state of an 
accreting gas and is of the order of unity, Cs is a sound speed of the gas and p 
is its density at infinity. 

When a particle is falling onto BH, some part of its energy is converted 
into radiation and the rest part increases the mass of BH. Hence the rate of 
accretion and the rate of BH mass growth is proportional to its luminosity. At 
the same time the accretion luminosity is limited by the Eddington condition — 
equality of the pressure of radiated photons and gravitational forces experienced 
by accreted plasma. This condition leads to the maximal value of luminosity 
for the stationary spherically symmetric accretion (the Eddington luminosity 
limit): 

. GAIm„c ,„qo M , 
Le = 4^ - 10^«— erg s~i, (4) 

(TT J-VIq 

where nip is the proton mass and ctt is the Thomson cross-section. 

Equation for BH mass growth in the case of stationary accretion of baryonic 
gas is 



dAI _ M 

dt tSed ' 

were the Salpeter time isai is 



(5) 



e Aic o e Le / v 

tsai = - — = 4.110«- ^ yr. 6 

1 — £ L 1 — £ L 

Here e is an efficiency of accretion and L is an observable luminosity of BH. It 
is also reasonable to suppose that L oc M as in the Eddington luminosity limit, 
see (|H). In this case the solution of equation (O is rather simple: 

M{t) = Afin exp[(< - <in)/<Sal] . (7) 
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where index "m" corresponds to the initial values. 

An accretion efficiency e is an emitted fraction of the total energy of accreted 
particle. To estimate this value, let us consider a particle with mass m orbiting 
around a nonrotating BH at some radius r. If the orbit is stationary, a total 
energy of particle is 

E^mc^ '- — . (8) 

{l-iGM/c^rf^ 

A difference between the energy at last stationary orbit ^ and those at infinity 
(r oo) is IS.E = 0.0572mc^. Respectively, for the extremely rotating Kerr BH 
this value is greater, /S.E = 0A2mc^. A real efficiency differs from an "absolute" 
value £ = 0.0572 for the Schwarzschild and e = 0.42 for the extreme Kerr case 
because some part of radiation is lost within the BH. In addition, some energy 
will be radiated during the particle motion between the horizon and the last 
stable orbit. For estimations it is usually used e ~ 0.1 for an accretion efficiency 
of nonrotating BH. 

A gravitational redshift z is the productive instrument in modern astro- 
physics. Let us denote a difference between the Newtonian gravitational poten- 
tials of two points as A$n- The observers located in these points will measure 
the different frequencies vi and V2 of the same electromagnetic wave. These 
values define the redshift: 

l + z= — ~1 + ^A$N. (9) 

An another group of formulae is related to the BH binaries and with the pro- 
cesses of their collisions and merging accompanied by an emission of gravita- 
tional waves. Investigation of gravitational wave generation in the Universe is 
the topical subject in view of coming completion of huge interferometric laser 
gravitational waves detectors. It is known that a collision of two BHs proceeds 
through a stage of quasi-periodic motion in close binary after their mutual grav- 
itational capture. An orbital period for two BHs with mass Mi and M2 is 



where a is a semi- major axis of the Keplerian orbit. Obviously, the gravitational 
waves emitted by binary have the same period. Corresponding duration of this 
Keplerian stage of gravitational waves emission by binary is [lOj 

5 a^c^ 

^"^^ " 256F{e) G3(Afi + Ma)^^' ^ ' 

Fie) = (1 - ey/' f 1 + + S^') > (12) 
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where e is an orbit eccentricity and ^ — MiM2/{Mi + AI2) 
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3 Observational signatures of SBHs 



Strictly speaking we have so far only numerous indirect indications in favor 
of BHs existence in the Universe. For this reason the claimed BHs of stellar 
mass in binaries and SBHs in galactic nuclei are only the candidates, not finally 
proven BHs. There is still some room for alternatives to BHs among suspected 
candidates. Decisive future experiment would be a demonstration or discovery 
of the BH event horizon, which is the most striking manifestation of strong grav- 
itational field and space-time curvature. At the present time, the crucial point 
of BH astrophysics is a revealing of observational signatures of event horizon in 
suspected BH candidates. The discovery of a BH event horizon would be the 
most crucial verification of the General Relativity as well. 

3.1 Active galactic nuclei and quasars 

Observations of host galaxies of quasars, i. e. the Quasi Stellar Objects (QSOs), 
are rather difficult because of the huge luminosity of central objects in compar- 
ison with stellar components and because of the large distances. Nevertheless, 
the connection between QSOs and their host galaxies has been very definitely 
established during the last decade though at least one exception is known - 
the bright quasar at the edge of a huge gas cloud without any visible stellar 
host galaxy llj. One of the explanation consists of a possibility that the SBH 
powering the QSO was ejected from the galaxy due to gravitational slingshot 
of three or more SBHs in the merger process of smaller galaxies [12] . The close 
similarity of the X-ray spectra of low- and highest-z QSOs means that SBHs 
with accretion disks were already formed z > 6. By using of some type of a 
"continuity equation" with "source" in redshift space it is possible to provide 
the link between SBHs in distant QSOs and SBHs in the centers of local galaxies 
[131 [21 [TSl [TB] . This analysis is suitable for a reconstruction of the accretion 
history of SBHs [T71[TH]. Different aspects of the accretion physics in QSOs are 
reviewd in [19] . 

3.2 Dynamics of gas and stars near SBH 

Specifical details of motion of gas and stars near the suspected SBH supply us 
with valuable information on the gravitational field in the vicinity of SBH. The 
best fit of numerous observations gives for a value of SBH mass at the Galactic 
Center Mbh ^ 3.6 • lO^Mo. 

Infra-Red (IR) imaging technic allowed to track the Keplerian orbits of 
several stars around the Galactic Center. Velocities of the stars approach 
9000 km s^^ (compare with the average velocity of stars 200 — 300 km s~^). 
Pericenter distances from the SBH are as small as 90 AU that is well inside a 
radius of gravitational influence of the SBH [20! . 

If a star passes close enough to SBH it could be ripped apart by the tidal 
gravitational forces. The conditions for such effect are favorable in the dense 
stellar nuclei of galaxies. The flare of UV- and X-radiation should be emitted 
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by the stream of stellar debris that plunges into the BH. These star destructions 
are really visible phenomena. Gezari et al. [21] observed the flare and its tail 
produced by stellar destruction near SBH in the galaxy at z = 0.37. Multi- 
wavelength data collected for 2 years beginning from the time of the event are 
in excellent agreement with theoretical predictions for the tidal disruption of a 
star. Gas streams from ripped stars serve the fuel for QSO activity. 

On the other hand, a vicinity of SBH is an appropriate area for the process of 
young star formation. Well known fact is the gravitational instability of massive 
accretion discs [35] . It could lead to formation of new stars if a cooling process in 
the disk is effective. Many young stars where discovered in the nearest vicinity 
of our Galactic Center, some of them about 0.03 pc from the center [53]. The 
processes of star formation and gas accretion could interfere meaningly [24] . 

In general, a final stage of matter accretion is poorly investigated yet because 
of inevitable observational difficulties in resolving small angular scales. An 
inner parsec of AGN contains dense molecular gas of H20, SiO and CH3OH, 
which is constantly heated by the accretion generated radiation. A detection of 
distribution of water maser emission line at 1.35 cm wavelength is a powerful 
investigation method of the near vicinity of SBH. A maser activity was first 
discovered in 1965, [25] . A water maser outside of our Galaxy was discovered in 
the spiral arms of the galaxy M33 by [55]. Up to now, in total about 60 events 
are found [57] . 

Observations of the broad fluorescence Ka iron line with energy 6.4 keV is 
also explained in the framework of accreting SBH in the distant galactic center 
[28] . Due to the combined Doppler effect in the fast rotating accretion disk 
and the gravitational field of the central BH the observed profile of emission 
line has an asymmetric double-peaked shape [29l [30l [311 (22] • The fitting of line 
profile corresponds to the emission of iron atoms at a distance of only a few 
gravitational radius of SBH. From the detailed spectroscopy of accretion disks 
it is possible also to evaluate the Kerr BH spin [33l [34l [35] [36] . 

New observational evidences continuously confirm the SBH paradigm. For 
example, it was reported recently [37] the activity of region close to SgrA* at a 
wavelength 3.5 mm. This region surrounds the central BH in the Galaxy having 
the size about 1 AU. On the other hand, the mass of the central object within 
this scale is well known, ~ 3 • IO^Mq. It permits to calculate an average density 
in the Galactic center, which appears to be ^ 6.5 lO^^Mo/pc'^ — too much to 
consists of gas in any state and/or stars. 

3.3 Galaxy — SBH correlations 

There are some kinds of observational data which must be explained by any 
model of SBH formation. On the other hand, these data give the "guiding 
thread" to choose the successful models. 

Early observations of normal galaxies, quasars and AGN recovered a simple 
relation A/sbh « L between a mass Msbh of the central SBHs in galactic nuclei 
and a total luminosity of host galaxies or their central stellar bulges [38l [Ml |4Q] . 
However an intrinsic data spread in the above linear relation was very high. 
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even larger then the observational errors. More tight correlations were found 
recently [41] between a mass Msbh and velocity dispersion erg at the bulge 
half-optical-radius: 

A^sbh«ct", (13) 

where a = 3.75±0.3. An other analysis 05] based on poorer set of observational 
data revealed the correlation with a somewhat different power index a = 4.80 ± 
0.54. In [43 it was shown that difference of slopes a arises mostly due to 
systematic difference in the velocity dispersions used by different authors for 
the same galaxies. More detailed analysis [331 of the set of 31 galaxies yields 
a = 4.02 ±0.32. Recently, a correlation between Msbh and host galaxy velocity 
dispersion for QSOs in the Data Release 3 of the Sloan Digital Sky Survey 
was discovered |44]. For small redshifts z < 0.5 these results agree with ones 
for the nearby galaxies. For a range 0.5 < z < 1.2 there is indication on the 
evolution with z: the bulges are too small for their central accreting SBHs. 
However, this evolution can be attributed to the observational biases. It have 
to be mentioned presumable correlation between the SBH mass and the age 
of the galactic stellar population [45]. This last type of correlation could be 
inconsistent with the scenarios of primordial SBHs origin. 

The observed correlations imply one-to-one relation between the processes of 
the central SBH and bulge formation and evolution in the all types of galaxies. 
Theoretical interpretation of these correlation is not very easy due to extremely 
different length-scales: the galactic bulge scale is a few kpc, whereas the scale 
of SBH gravitational influence is millions times smaller. Specific mechanism 
is needed for mass transfer from the bulge to its innermost part. One of the 
possibilities is that SBH formation depends on the dynamical properties of bulge 
through the rate of gas supply. In this case a velocity dispersion a measures the 
depth of the gravitational potential in which a SBH forms and grows. Another 
discussed possibility is the energy feedback between an accreting SBH and bulge 
which results in a self-regulated accretion and a corresponding growth of the 
central SBH [JB] . This model predicts the following relation, Mbh oc cr^ ■ Third 
model [17] supposes that bulges and SBHs were formed during the multiple 
galactic merging. The resulting correlation between Mbh and CTe in the galactic 
merging model is in a good agreement with observations. A forth model proceeds 
from the hypothesis that observed correlations are stochastic in origin [48] . 
The basic assumption of stochastic model is an existence of a pre-galactic or 
primordial SBH population. A similar hypothesis of an existence of pre-galactic 
SBH population was used [HI I208j for interpretation of the QSO evolution. 

3.4 Black holes of intermediate mass 

Some mechanisms of SBH formation imply the Intermediate Mass BHs (IMBHs) . 
The latter are believed to be an intermediate stage of SBHs growth in astro- 
physical evolution scenarios with an accretion and merging. Fast variability of 
compact X-ray source in the starburst galaxy M82 indicates that this source is 
an accreting BH ^50j. With the observed luminosity of lO'*^ erg s^^, a mass of 
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the supposed BH is of the order of ~ lO^M© provided the emission is isotropic 
and corresponds to the Eddington Umit. The IMBH in M82 exists in star burst 
environment. This is a hint on a possible formation of the IMBH in the M82 by 
a runaway merging of stelar mass BHs ^T\. An other suspected IMBH in the 
star cluster Gl of the Andromeda galaxy has mass 2 • lO^M© [52], but there is 
still a room for a tight cluster of massive stars in this case. 

Some of nearby Globular Clusters may contain IMBHs according to the 
measurements of stellar kinematics. On the basis of detailed observations of 
Globular Clusters it was stated that IMBHs may inhabit about a half of Globular 
Clusters in the Galaxy [53]. These Globular Clusters are very old in general. 
In particular, the mass of IMBH containing in M15 is about 4OOOM0 [M]. It 
is interesting that masses of suspected IMBHs in Globular Clusters are almost 
those as expected from the extrapolated AfgeH-c correlation for galaxies, in 
spite of poor statistics of the central BHs with AIsbh < lO^M© in galactic 
nuclei. 

It was shown [5^ that the primordial IMBHs could be connected with such 
dense Dark Matter (DM) objects like neutralino stars through their common 
origin from primordial density perturbations. The neutralino stars with mass 
~ (0.01 1)M0 were proposed to account for observed microlensing events in 
the Galactic halo [SS] . These neutralino stars consist of weakly interacting DM 
particles, and were formed in the expanding Universe from adiabatic density 
perturbations. The existence of neutralino stars should give rise to a large num- 
ber of primordial IMBHs with a mass ~ IO^Mq. These IMBHs may constitute 
a large fraction of the Galactic halo mass. 

3.5 Prospects for gravitational waves detection 

Coalescence of BHs in the galaxies and clusters are inevitably accompanied by 
strong bursts of gravitational radiation. The planned ES A/NASA space obser- 
vatory LISA will be capable to detect these coalescence events in very distant 
galaxies. The combined theories of galaxies and SBH merges are developed to 
predict the principally observed gravitational waves signals [57] jSH] [59] . Future 
observations of binary merging with LISA will allow to test General Relativity 
in details, place bounds on alternative theories of gravity and study the merger 
history of massive BHs [50] . 

3.6 Astrophysical constraints on SBHs 

There are several purely astrophysical constraints on the mass and number of 
SBHs. If SBHs provide the dominant part of DM in the Galaxy, then they must 
tidally interact and disrupt some part of Globular Clusters. The SBH mass was 
constrained in this case, Mbh < lO'^M© [61]. At fieH ~ 1, primordial BHs are 
capable of distorting the CMB spectrum if they are formed about 1 s after the 
annihilation of e"'"e~-pairs [55]. Mass accretion at the pre-galactic and present 
epochs contributes to the background radiation in different wavelength ranges. 
However, calculations are strongly model dependent and yield Hbh < 10~^ 



9 



lO"'^ for A/bh ~ IO^Mq. In (53], the constraint on the fraction of intergalactic 
BHs, r^BH < 0.1, was obtained from a condition of absence of the rehable 
gamma-ray-burst lensing events for lO^M© < Mbh < lO^M©. A more stringent 
lensing constraint, f^BH < 0.01 for the mass range lO^M© < Mbh < lO^M©, 
was obtained from VLBI observations of compact radio sources |M] . 

3.7 Dark matter spike in the Galactic center 

Annihilation rate of weakly interacting cold DM particles at the Galactic Center 
could be greatly enhanced by the growth of density spike around the central 
SBH. The resulting annihilation signal may be observed in the form of high- 
energy gamma-rays, positrons, antiprotons and radio emission |65| . It could 
facilitate the study of DM content. At the same time this annihilation spike 
may be the observational signature of the SBH in the Galactic center. The 
form and amplitude of this spike is rather uncertain due to influence of different 
dynamical processes such as sinking of baryonic gas and hierarchical merging of 
sub-halos |661 167j . Searching for annihilation signals could set an upper bound 
to the DM cusp, DM particle annihilation cross-section and the merger history 
of the Galaxy. The DM spikes could also arise around primordial IMBHs |3S], 
and the DM annihilation in these spikes is in observable range of future detectors 
[68] . A model of adiabatic growth of a seed BH in coUisionless DM halo through 
the DM accretion is proposed in ^55] . The supposed combined evolution of SBH 
and DM halo results in the observed Msbh — cr relation. 

4 Astrophysical scenarios of SBH formation 

4.1 Dynamical evolution of galactic nuclei 

Physically the most simple and natural possibility for SBH formation is a dis- 
sipative dynamical evolution of a dense stellar cluster in the galactic nucleus 
ending by its total gravitational collapse. This idea was first proposed at the 
dawn of quasar discovery and then was elaborated in some details in numerous 
works. 

An evolved galactic nucleus contains dense central stellar cluster and prob- 
ably subsystem of compact stellar remnants, such as neutron stars and stellar 
mass BHs. This subsystem come from the dynamical evolution of star cluster in 
the galactic nucleus through merging of stars, thereby forming the short-living 
massive stars. The clusters of compact remnants form almost inevitably. 

To perform estimations, let us consider a central stellar cluster in the galactic 
nucleus of mass M and radius R (the radius can be effectively defined through 
the sphere which contain half of the cluster mass), consisting of TV ^ 1 of 
identical stars with mass m = M/N and radius r* . The virial theorem gives the 
velocity dispersion in the cluster v ~ {GM /2RY/^ . 

A characteristic measure of the stellar system compactness is a depth of 
its gravitational potential |^| ~ . That is why the value of star velocity 



10 



dispersion v is crucial parameter in evolution history of the stellar system. In 
the case v c the Newtonian approximation is applicable to a star motion in 
cluster. The Newtonian approximation is failed when a stellar system becomes 
relativistically compact, with |(/)| ~ c^. At this stage of evolution the radius of 
the systems is close to its gravitational radius. The cluster as a whole or its 
major part becomes dynamically unstable and collapses into the SBH. 

Dynamical evolution of the cluster at nonrelativistic stage was widely dis- 
cussed (see e. g. [TDl [ZIl HH IZl] and references therein) and almost exhaustively 
studied in the Fokker-Plank approximation. It turns out that very important 
factors of the evolution are the tidal interactions of stars, the binary stars for- 
mation and the star mass segregation. Due to the mass segregation, the binaries 
and most massive stars are concentrated in a central part of the cluster. 

Dynamical evolution of stellar systems those as galactic nuclei which are 
dense enough, proceeds through the stellar collision and coalescence stage [73l 
[111 [75l |72l [76] . This process is enhanced by the formation and hardening of 
star binaries formed by tidal two-body interactions [771 [ZSl [75] • The newly 
formed hard binaries result in heating of the cluster. The resulting evolution of 
star cluster is determined by the competition of tidal dissipation processes and 
heating by newly formes binaries. In the clusters with a high enough velocity 
dispersion v dissipative processes dominate and accelerate the contraction and 
growth of V. Later, due to the growth of dispersion v the stellar collisions 
become catastrophic and destroy the colliding stars. As a result the compact 
stellar cluster transforms into the Supermassive Star (SMS) with an appreciable 
fraction of neutron stars and and BHs of steUar mass [80 l [ST ] [82 l [83 l [M l [85] . 

There are two distinct ways of IMBH or SBH formation in the evolving star 
cluster: 

(i) Central core of the cluster reaches relativistic state and collapses as a 
whole. 

(ii) Stellar mass BHs and neutron stars merge (aggregate), progressively 
increasing their mean mass. 

In the following we suppose for simplicity that stellar cluster consists of 
identical stars of mass m. The main virial parameters of a self-gravitating 
stellar system are a total mass M — Nm, a mean radius R within which N/2 
stars are contained, a star velocity dispersion is v. These parameters obey a 
virial theorem |86| : 

GmN 



,.2 



V 



2R 

A total energy of the self-gravitating system is 

1 



(14) 



E = --mNv\ (15) 

These virial relations are established in the dynamical time idyn = R/v. The 
corresponding dynamical evolution equation on the time-scales exceeding tdyn 
is derived by the differentiation of viral relations ([T4l) and ([TS]) : 



R E 

— = h2— . (16) 

REM ^ ' 
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A central stellar cluster in a galactic nucleus must be compact enough to evolve 
during the life-time of the galaxy. Due to the Coulomb character of the gravi- 
tational interactions of stars the evolution of star orbits in the cluster proceed 
diffusively with a characteristic two-body relaxation time [571 IHH] 

2X 1/2 



S^GWjiA' ^^^^ 

where m is a constituent stellar mass, n is a local star number density, A ~ 
log(0.4A^) is a gravitational Coulomb logarithm. For many-body stellar systems 
with number of stars ^ 1, the relaxation time is much greater than a 
dynamical time tdyn- 

— ~ r^TT?: > 1 at A^ > 1. (18) 

tdyn log(A^) ^ ^ 

This means that stars in the self-gravitating dynamically stationary systems 
have regular orbits with integrals of motion (energy and angular momentum) 
slowly changing at time-scales t <^t^. 

The smallness of the ratio of relaxation time and the age of the galaxy 
Tg ^ is a measure of stellar cluster compactness needed for a substan- 

tial evolution. Stellar clusters with t^/Tg ~ <^ 1 are called (dynamically) 
evolved: they have enough time for substantial change of initial orbital energies 
and angular momenta of stars in the cluster due to two-body interactions. The 
examples of dynamically non-evolved stellar clusters with t^/Tg ^ H^^ ^ 1: 
stellar disks of spiral galaxies. Globular Clusters, galactic DM halos. The con- 
ditions for substantial evolution are realized only in the central most compact 
parts of these objects. For the same reason the central parts of the Globular 
Clusters and galactic nucleus are suspected for the final gravitational collapse. 



4.1.1 Non-dissipative contraction of star cluster 

A non-dissipative dynamical evolution is governed by the evaporation of fast 
stars from the system. The rate of fast star evaporation due to two-body inter- 
actions is 

A^ 

A^ = -aev — , (19) 

where numerical coefficient a^v — 10^^ in the Fokker-Plank approximation [701 
I71j . Due to diffusive nature of star orbital energy changing, an evaporating 
star escapes from the stellar system with a nearly zero total energy, i.e. with a 
nearly parabolic velocity v^sc = 2u. Hence the total energy of the stellar system 
is nearly conserved during star evaporation, E ~ const. This defines the scaling 
of the virial parameters: the evolutionary growing of a star velocity dispersion 
with a diminishing of the number of stars, v cx N~^^^, and a corresponding 
diminishing of the stellar system radius, R (x A^~^/^. Integration of evolution 
equation ([T6|) at the condition E = const by using p9|) and by neglecting 
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variation of logarithmic factor A gives a simple law for non-dissipative dynamical 
evolution 



where the lifetime of the system t^v — (2/7acv)ir(0) and tr{0) is a relaxation 
time at the initial moment t — 0. According to this formula, a contracting 
stellar system reaches the relativistic stage, v ^ c. At the relativistic stage, 
not the evaporation of stars but the dissipative processes such as gravitational 
radiation and star collisions become the crucial dynamical evolution factor. The 
dissipation processes hasten the contraction of cluster towards the onset of dy- 
namical instability and final gravitational collapse. This scenario of evolution 
is valid in clusters consisting of degenerate stars like neutron stars and/or BHs 
of stellar masses. In the case of contracting cluster of normal stars, the dissipa- 
tive tidal interactions and the star collisions become more important than star 
evaporations long before the relativistic stage of evolution. 

4.1.2 Dissipative contraction of star cluster 

The influence of dissipative tidal interactions and star collisions on the dynam- 
ical evolution is crucial in the star clusters with a large value of stellar veloc- 
ity dispersion v, namely when v ^ Vp, where Vp — {2Gm/r^y/^ is a surface 
parabolic velocity of constituent stars. For example, in the case of the Sun 
Vp ~ 620 km s"^ < c. 

In close encounters the non-radial oscillations are excited by tidal interac- 
tions of stars. Some part of kinetic orbital energy transferred into these oscil- 
lations. Finally the energy of non-radial oscillations is thermalized inside stars 
and reradiated. This is also the mechanism for dissipative tidal friction of stars 
and the tidal capture mechanism of close binary star formation |77j . The orbital 
energy deposited into oscillations during one close encounter equals [79j 



for a specific case of stars with n — 3 polytropic structure. Here i?min is a 
periastron distance. The resulting tidal "cooling" rate of the stellar system in 
the case of the Maxwellian distribution of stars over velocities is |Sni HB US 



In this equation r(0.9) ~ 1.069 is the Gamma-function and — is a 

mean kinetic energy of stars in a cluster (or cluster's "temperature"). The tidal 
dissipative friction of stars in the cluster accelerates its contraction. Under the 
influence of both evaporation of stars and tidal friction the evolution equation 
p6|) according to (|19p and ((22)) takes the form 




(20) 




(21) 




(22) 
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where 
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3(V2)V^aevA 

r(o.9) 



-, 5/9 



Vdis = 



•p- 



(24) 



For a case of cluster consisting of the Sun type stars, Vdis — 190 km s^^. Ac- 
cording to equation (I23p dynamical evolution follows non-dissipative evaporative 
scenario ([20| until v <^ Wdis- It also concerns the central parts of the nowadays 
Globular clusters with v ^ 10 km s^^. Only a very small fraction of stars in 
Globular Clusters survives in the final steps of dynamical evolution. Addition- 
ally a lot of hard binaries (with binding energies exceeding are formed 
in Globular Clusters by dissipative two-body encounters ^7\. As a result the 
heating of cluster by these hard binaries [78j prevents the formation of a massive 
BHs in the Globular Clusters (see for details 

On the contrary, in stellar systems such as the galactic nuclei with v > vns, 
the tidal friction dominates over star evaporation and heating by hard binaries. 
As a result the cluster contracts formally to a zero radius having a finite number 
of stars in the remnant [90l |92] . Physically this dissipative stage of evolution is 
finished by the direct collisions, merging and destruction of stars. Collisions of 
stars with relative velocity u < Vp result in merging of stars or their explosion of 
a supernova type. A head-on collisions of stars with u > Vp are more dangerous 
and inevitably end with a total destruction of these stars 73 . In this limit the 
rate of disruptive collisions of stars in the system is equal to 



A contracting star system with v > Vp consists of a decreasing number of stars 
and an increasing number of moving gas clouds — remnants of destructive 
star collisions. These gas clouds participate in the virial balance as the peer 
entities until they are small compared with the size of the remaining system. 
In a simple approximation, suppose that the two clouds with a nearly equal 
masses are formed as a result of the two stars catastrophic collision, and so 
Nm — M — const. Integration of equations (|16p and (|25p in this approximation 
gives 



where the collision evolution time icon — [vp / v{0)]'^ At,-{Q) . The system contracts 
according to ([^ until its radius reaches R ~ N^/^r^,, at which point tcoii ~ 
tdyn- At this point the system ceases to exist as a collection of separate stars 
and transforms into the compact self-gravitating massive cloud of gas — the 
Supermassive Star (SMS) [9l El [95l |98]. The possible evolution tracks of 
stellar clusters are shown in the Fig. [T] 

4.2 Collapse of supermassive star 

A Supermassive Star (SMS) was first proposed long ago as a possible source 
of quasar activity and nowadays is considered as a short intermediate stage of 




(25) 




(26) 
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Figure 1: Evolution tracks of stellar systems in the plane "number of stars N, 
velocity dispersion w" (see [92] for details). A — globular clusters, B — normal 
galaxies with low v, C — galactic nuclei with high v In region I, the tidal 
dissipative effects and heating of stellar system by newly produced binaries do 
not influence the evolution of the system which follows to the non-dissipative 
evaporative scenario (see Section [4 .1.11) . The boundary between regions I and 
II (and between III and IV) corresponds to condition E = 0, when dissipative 
effects in the system are compensated by the heating by binaries. In the region 
II, the heating of the system by newly formed binaries dominates over the tidal 
dissipation, E > 0. The boundary between regions II and III is determined 
by the condition v — 0, or equivalently v — t^^J^. In the region III, where 
V > ~ 190 km s~^, the heating by binaries prevails over the tidal dissipation 
and V decreases with the decrease of N. In the region IV, the evolution of 
stellar system is mainly affected by tidal dissipative energy losses which lead to 
a rapid contraction of the system toward to the final collapse into SBH. Above 
the dotted line the rate of star collisions becomes greater than the rate of star 
evaporation. The dashed-dotted line corresponds to — 10^° yr. 
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galactic nucleus evolution toward the formation of SBH. To prevent a fragmen- 
tation before a SMS formation, a primordial gas must be hot and magnetized. 
[HI [Ml 1100] . The natural way for SMS production is stellar collisions stage of 
evolution of compact galactic nuclei [TH [HUl H] ^ . 

Fast evolving SMSs include an unstable hydrodynamic relativistic radial 
mode [SSI IfOll 11021 11031 1104] and eventually collapse to form one or several 
close SBHs [94j [95l [98l il05| . Recently numerical simulations provides the possi- 
bility to analyze the final state of the SMS gravitational collapse [1061 11071 1108] . 

The SMSs may be naturally formed in the galactic nuclei from gas produced 
in the destructive collisions of stars in the evolved stellar clusters with a velocity 
dispersion v > Vp. Characteristic time-scale for the dynamical evolution of 
stellar cluster is the (two-body) relaxation time tr- According to ([17]) it can be 
expressed as 

tr ~ 4.6 X lO^A^I (j^^ yr, (27) 

where N — lO^A^g is the number of stars in the cluster. The corresponding virial 
radius of the stellar cluster is i? = GNm/2v'^ . At v > Vp, where Vp is an escape 
velocity from the surface of constituent star, the time-scale for self-destruction of 
normal stars in mutual collisions is tcoU = {vp/v)'^Atr [72] . Numerical modelling 
of catastrophic stellar collisions has been performed by e. g. [1091 1110[ . We 
choose w ~ Up as a characteristic threshold value for a complete destruction of 
two stars and final production of unbound gas cloud. As a result, the stellar 
cluster in the evolved galactic nucleus with v > Vp transforms into the SMS due 
to catastrophic stellar collisions. At v :^ Vp the time-scale for the formation of 
SMS due to destructive collisions of stars is of the same order as the relaxation 
time, icoii(wp) ^ ir(wp). A total mass of the gas produced by destruction of 
normal stars composes the major part of a progenitor central stellar cluster 
in the galactic nucleus. The natural range of masses for the formed SMS is 
MsMS = 10^- IO^Mq. 

A newly formed SMS with mass Msms and radius i?sMS gradually contracts 
due to radiation with the Eddington luminosity, L = Le- A nonrotating SMS is 
a short-lived object that collapses due to post-Newtonian instability. Rotation 
provides the stabilization of SMS if the rotation energy is an appreciable part of 
its total energy i?sMS — {GMgj^g/2RsMs)- In general an evolution time of SMS 
is the Kelvin-Helmholtz time-scale ^sms = ^'SmsZ-^e oc Rg-^g- Stabilization 
of SMS by rotation (and additionally by internal magnetic field) ensures in 
principle its gradual contraction up to the gravitational radius [98l 11111 1112| . 
A resulting maximum evolution time of SMS is of the order of the Eddington 
time = O.lMsMSC^/^E — 4.5 x 10'' yr. The distribution of gas in a SMS can 
be approximated by the polytropic model with an adiabatic index 7 — 4/3. For 
this value of adiabatic index the central gas density in SMS is pc = fcc'T-SMS'Ttp 
and central sound velocity Cs,c = ^s^smS: where nsMS is a SMS mean number 
density, usms = (GMsMs/2i?SMs)^^^ is a SMS virial velocity and numerical 
constants kc ~ 54.2 and ks — 1.51 respectively. 

A SMS, formed in the galactic nucleus, may contain compact remnants of 
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exhausted stars in the form of neutron stars and stellar mass BHs. Due to 
dynamical friction these compact objects are settled down to the central part of 
the SMS. They form a very compact and fast evolving self-gravitating subsystem 
which collapses into the massive BH earlier than the host SMS. An observational 
signature of these nearly collapsing SMS in the distant galactic nuclei may be 
a long-wavelength gravitational radiation |113j and a power flux of high-energy 
neutrino [841185] . 

4.3 Collapse of neutron star cluster 

Neutron star cluster evolves much faster than a host SMS and collapses finally 
into the massive BH. The dense clusters of neutron stars and stellar-mass BHs 
can undergo the catastrophic events of gravitational collapse due to general 
relativistic instability. The collapse goes through the avalanche-type contraction 
described firstly in [114] and confirmed by the detailed numerical calculations 
of [m ia IMl Hill HIII [TTTI inij . They found that even in the case of extremely 
centrally condensed cluster configuration with an extensive Newtonian halo, 
a significant fraction (several percents) of total mass collapses into a central 
BH in a few dynamical times. The collapse begins when the central redshift 
approaches 2 ~ 0.5 and a corresponding virial radius of the cluster diminishes 
up to i? ^ 3Rg. The 'avalanche' arises because the falling of stars onto the 
center leads to an increase of the gravitational field. The latter acts on the other 
particles/stars, whose orbits contract in turn, and so on. The non-circularity of 
star orbits are of principle significance because they connect the different layers 
of the cluster. Only stars at highly elongated orbits are involved in the collapse 
and accreted into the growing BH. After the collapse the cluster settle down to 
a new, dynamically stable equilibrium state: central BH embedded into quasi- 
Newtonian cluster of stars and compact stellar remnants. This evolutionary 
track is time-consuming because it requires the prior evolution of normal stars 
and subsequent settle down of their remnants to a highly concentrated system. 

The further growth of a seed massive BH with a mass A/bh > W^Mq is 
governed by an accretion of ambient stars, gas and DM particles (see e. g. 
pi ll20l[l21| V During a lifetime of the normal galaxy a seed massive BH should 
grow sig nificantly, up to Mbh ^ 10^ - IO^Mq. 

4.4 Merging of stellar-mass black holes 

The alternative model of IMBH formation in a star cluster is the multiple merg- 
ing of stellar mass BHs when. The evolution of mass distribution usually studied 
by using the Smoluchowski equation fSH I122[ 1123] . Under the conditions that 
may exist in some star clusters the characteristic timescale for runaway merging 
is ~ 10^ yr, which is short enough for the formation of BH with mass ^ 10'^Mq. 
More massive BHs sink into the cluster center under the influence of dynamical 
friction and increase the merger probability in the core. Such a 'mass segre- 
gation' is very important process for the evolution of the cluster core. After 
merging of hosted protogalaxies, the IMBHs may sink to the galactic centers 



17 



and merge into SBHs. It should be mentioned that this evolutionary track is 
rather time-consuming and probably can not explain the observed early QSO 
activity. Additional prolific accretion of baryons or DM is required [124j . It is 
considered also a possibility of the super-Eddington accretion regime, [1251 1126j 
which strongly accelerates a BH growth. 

A peculiar example is an accretion onto BHs of the Dark Energy (DE) with 
an equation of state P{p) < 0. In the extreme case of phantom energy (with 
P + p < 0) , the accretion is accompanied with a gradual decrease (rather than 
increase, as in a case of the usual matter) of the BH mass |127[ 1128] . Re- 
spectively, masses of all BHs tend to zero while the phantom energy Universe 
approaching to the Big Rip. 

4.5 Merging trees of galaxies and SBHs 

Galaxies are formed by hierarchical merging of smaller protogalaxies, and even 
nowadays a significant fraction of low redshift galaxies still experience merging. 
When galaxies merge, their cental SBHs get a chance to merge too. The dynam- 
ical friction in the field of stars and DM brings the SBHs into the central part of 
coalesced galaxies, where they finally merged into a bigger SBH. The timescale 
for a 1O®M0 SBH to spiral down into the nucleus of typical giant galaxy from 
an initial radius of 10 kpc is ~ 10^° yr. A subsequent closing of two SBHs to 
one another proceeds even faster ^ 10^ yr due to a high number density of stars 
in the galactic nucleus, which produces a strong dynamical friction ^1^51 1131] . 
On the final stage of merging the gravitational radiation dominates. However, 
there are calculations which predict the long living binary SBHs, longer than 
the Hubble time [132 . 

Merging of small protogalaxies containing the IMBHs is a promising mecha- 
nism of SBH formation. The seed IMBH could be formed initially in the mini- 
haloes collapsing at z ~ 20 from high-cr density fluctuations [133|, I134|, I135j . 
They increase their masses up to modern values 10^ — IO^Mq due to multiple 
merging and accretion |136[ I137j . Direct evidence of SBHs merging was founded 
by VLA telescope from observation of jets flipping (X-shape of radio lobes) in 
NGC 326 galaxy [143j . It is interesting that a total energy transferred from the 
binary SBH to the field stars is comparable with a binding energy of galactic 
core [129i . A combined upper limit on the mutual evolution of SBHs and their 
host galaxies was obtained in |130j . 

Recent X-ray observations clearly demonstrate the existence of binary AGN 
in the galactic system Arp 299 consisting of two galaxies in an advanced merging 
state, NGC 3690 and IC 694 to the east, plus a small compact galaxy |138j . The 
other cases of two AGN in a merging system is NGC 6240 [l39] and NGC 6104 
|140j . A close encounter and coalescence of spiral galaxies triggers an infiow of 
nuclear gas, which fuels both a powerful starburst and strong accretion onto the 
central SBH [l4n[T42] . 

Coalescence of SBHs in galaxies must be accompanied by the strong bursts 
of gravitational radiation. The planned ES A/NASA space observatory LISA 
will be capable to detect these coalescences and verify model predictions for 
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SBHs merging rate [1431 1144] . A binary SBH leaves an imprint on a galactic 
nucleus in the form of a "mass deficit" , a decrease in the mass of the nucleus 
due to ejection of stars by the binary [145j . Comparison with observed mass 
deficits implies between 1 and 3 mergers for most galaxies, in accordance with 
the hierarchical galaxy formation models. 

A mean number of galactic mergers is rather low to ensure the growth of 
stellar mass BHs to SBHs in most galaxies with the only exception of giant CD 
galaxies in the centers of rich clusters. Another difficulty is the "final parsec 
problem": quite large distance between SBHs in binaries after the merging of 
host galaxies to provide the strong enough gravitational wave radiation losses for 
merging of SBHs in binaries [146] . Nevertheless, it seems that SBHs populate 
near the all spiral and elliptical galaxies. This may be a hint in favor of some 
universal mechanism of SBH formation long before the formation of galaxies 
in the early Universe. The merging rate of galaxies and their central SBHs is 
insufficient to ensure the needed mass growth at least at low redshifts, z < 0.36 



5 Cosmological scenarios of SBH formation 

The present state of art in understanding of SBH origin provides only a list 
of different scenarios, more or less elaborated. Even Moreover, we even non 
confident in the choosing between two major alternatives: the galactic or the 
cosmological origin of these objects. In view of these uncertainties we will 
itemize in this review the most popular known scenarios with a brief elucidation 
of underlying physics. It is a challenge for future theoretical investigations and 
detailed astrophysical observations to recover the real mechanism(s) of SBH 
formation. 

The problem of SBH origin becomes more acute with a discovery of very 
distant quasars. Most distant of them are: 



Observations of high redshift quasars lead to an unambiguous conclusion that 
host galaxies of SBHs are formed very early. The existence of quasar with 
redshift z = 6.43 indicates that SBHs with mass around 1O^M0 already existed 
when the Universe was less than 1 Gyr old. In order of magnitude, it is just 
a lifetime of an ordinary massive star. Such stars produce BHs with masses of 
order IM0 after explosion. Using formula ([7]) one can obtain the time needed a 
BH to grow up in 10^ times. In case of the Eddington accretion, a BH increases 
its mass in 10^ times during 800 million years. There is no room for first ordinary 
star formation. 

Nowadays we have a list of competing models to solve this puzzle. This 
list consists of two parts: cosmological (pre-galactic) astrophysical (galactic) 



J114816.64+ 525150.3, 
J103027.10 + 052455.0, 
J130608.26 + 035626.3, 



z = 6.43,M~31O^M0; 
z = 6.28, Af~ 3.6 IO^Mq; 
z = 5.99,Af ~ 2.41O^M0. 



(28) 
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models. The latter were discussed in previous Section. Before the discussion 
of cosmological models let us consider a couple of "intermediate" models where 
SBHs and galaxies are formed at the same time approximately. 

5.1 Population III stars or supermassive stars 

BHs with masses about IOOMq are thought to be produced as remnants of first 
stars (Population III stars) which are supposed as massive as 2OOM0 |148|. I149|, 
11501 1151] . A lifetime of these stars is quite short — after several millions years 
they collapse into BHs. 

According to a hierarchical model of structure formation, less massive DM 
objects are formed earlier. Formation of gravitationally bound objects with 
^ 10^ — 1O^M0 is expected to be before z ^ 10. The short-lived massive stars 
with mass ~ lOOM© can appear from the metal-free gas clouds in the first DM 
sub-haloes, with a several massive stars per a DM sub-halo. 

Formation of population III stars requires an effective cooling mechanisms 
at early epoch of fragmentation of primeval gas clouds into the stars. Usually 
believed that the line emission of molecular hydrogen is the basic cooler in these 
primeval gas clouds. If the gas cloud evolves without fragmentation to stars, it 
can reach the state of a single Super-Massive Star (SMS) with mass in the range 
^ 10^ — IO^Mq. a SMS is a very short-lived object and inevitably collapses 
into an IMBH. 

5.2 Collapse of massive primordial clouds 

SBHs could be formed as a result of dissipation and collapse of primordial gas 
during the early stages of galaxy formation [6l[T4 t[T52lll63|.[T64jll65|.fT66] . Initial 
masses of such SBHs are of order IO^Mq. Another scenario of SBH formation 
in the low angular momentum haloes was proposed and recently considered in 
details by [152j . In this scenario the baryons cool and sink to the center of the 
halo and finally settle down into a rotating disk. Due to viscosity, an angular 
momentum is transferred outward from an inner region of the disk. A model 
of the formation of lightweight BHs (a few solar mass) by the direct collapse 
of gaseous clouds (without a preceding formation of stars) was elaborated in 
|153j . It is supposed that these BHs are growing up to the supermassive ones 
by accretion of massive envelops with the super-Eddington rate. 

The haloes are distributed over their angular momenta. This distribution 
depends only on the spectrum of primordial density perturbations and was ob- 
tained from analytical calculations and numerical simulations [1541 11551 11561 
11571 1158] . Only those haloes from low angular momentum tail of the distribu- 
tion are capable to collapse. According to [152j. haloes more massive than a 
critical threshold about ~ 7- IO^Mq at red shifts z ~ 15 contain the gravitation- 
ally unstable discs with an efficient viscosity. The authors of [152] also showed 
that under reasonable assumptions, this model leads to the observed correlation 
between BH masses and spheroid properties. 
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Combined dynamics of baryons and DM in the halo of forming galaxy was 
analyzed in [159| . The cooHng and contraction of baryonic component both lead 
to the formation of SBH which grows rapidly due to accretion as baryons [160j 
and the DM [TOTHTB^ . Absorption of DM in [TBT] and [TO^ was examined by 
taking into account the scattering of DM particles on stars near the galactic 
center. It was shown that significant flux of DM on a seed BH in the galactic 
center of galaxy dominates. Respectively, a fraction of DM in the total mass of 
the SBH in this scenario is significant. 

5.3 Pre-galactic SBH 

As it was mentioned above there were proposed some mechanisms of massive 
BH production before the formation of first stars. Below we shortly list the 
proposed mechanisms of primordial BH formation (see also the detailed review 
|167j ) and discuss some of the listed models in more details in the next sections. 

Adiabatic fluctuations at radiation-domination stage. According to 
[62l [96l \97\ I168j , density fluctuations at the radiation-dominated evolutionary 
stage of the Universe give rise to primordial BHs. The formation of BHs is 
quite efficient in the case of "blue" spectra of power-law density perturbations 
n > 1.1 - 1.2. 

Isocurvature (or isothermal) fluctuations after inflationary stage. 

Variety of models with the flat or bumped spectra of isocurvature fluctuations 
were proposed [Hi HTDl [IZIl HZl HZl HZl ^JE HZi \I77\ A very promising 
mechanism of SMB formation from the large amplitude isothermal fluctuations 
in baryonic charge density was proposed by Dolgov and Silk |178j . In this 
mechanism the bumped spectrum of isothermal fluctuations is a byproduct of 
vacuum bubble collapses during phase transition at the inflationary stage. 

Artificial form of the inflaton potential. Models of such sort represent 
a modification of the previous model. It is known that the inflationary stage is 
the reason of density fluctuations. Spectrum of the fluctuations depends on the 
form of the inflaton potential, see i. g. jl79[ 1180] . As it was shown in fl81', it 
opens a possibility to produce high density fluctuations that collapse into SBHs 
in the following. 

SBHs ft-om phase transitions. BHs with masses ~ IMq possibly formed 
at quark- hadron phase transition at the cosmological moment 10~^ s, (18211183] . 
Such BHs would be a component of DM today. 

Soft equation of state. Suppose that some massive non-relativistic par- 
ticles dominates at an early period of the Universe evolution. In that case the 
pressure is negligible and could not prevent gravitational collapse of high density 
regions [168]. 

Bubble collision. False vacuum decay is usually accompanied by formation 
of spherical walls with a true vacuum inside them [IM IT851 [TM ITSTl fT88l IT891 
11901 1191j . The walls quickly expand and collide what could lead to collapse of 
islands of false vacuum into BH. Most massive BHs of order IMq are produced 
during the period of quark-hadron phase transition. 
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Closed domain walls. A mechanism of BH formation from the closed 
domain walls was proposed in [192[ 11931 1194j . These domain walls could be 
originated due to evolution of a scalar field during inflation. An initial non- 
equilibrium distribution of scalar field imposed by the background de-Sitter 
fluctuations gives rise to the spectrum of BHs, which covers a wide range of 
mass — from the subsolar up to supermassive ones. The primordial BHs of 
smaller masses are concentrated around the most massive ones within a fractal- 
like cluster. It was revealed that this mechanism is a rather common for many 
inflationary models and it worth to discuss its essential details at the end of the 
review. 



5.4 Adiabatic fluctuations at radiation-domination stage 

Gravitational collapse and the formation of a primordial BH occurs if a relative 
radiation density fluctuation Sn = {p — Pc) / Pc satisfies some conditions |62i I96i 
[97] at the moment t of its separation from the cosmological expansion. Namely, 
the BH is nucleated if 

(Jc < < 1, (29) 

where 5c — 1/3 [62]. The left-hand inequality implies that the radius of the 
perturbed region at the time t exceeds the Jeans radius ct/\/i. The right- 
hand inequality corresponds to the formation of a primordial BH rather than 
an isolated Universe. The primordial BHs are produced with a near horizon 
size, and thus their masses are connected with the formation time as 

GM f M \ 
t = ^- = 0.5 — s. (30) 

In recent years, numerical experiments have revealed a near-critical gravitational 
collapse of primordial perturbations with less than the horizon scales |1951 11961 
11971 1198] . In this process a mass of forming primordial BH is 

AfBH = AAfH(fe-^c)^ (31) 

where ^ ~ 3, 7 ~ 0.36, and 5c — (0.65 0.7). The mass (|3T|) may be much 
smaller than the horizon mass Afn. However, as shown in ^199] . the primordial 
BH mass distribution for critical gravitational collapse is concentrated near 
AfsH ~ Afn, and the contribution of low masses to the cosmological primordial 
BH density is modest. 

If the power spectrum of primordial cosmological perturbations is a power 
law with an index n > 1, then the primordial BHs are formed in a wide range 
of masses. However a power law spectrum with n > 1 is not favored by recent 
CMB observations, and a corresponding number density of resulting primordial 
BHs is extremely (exponentially!) small. This is because the primordial BHs 
form at the tail of the Gaussian distribution with probability 

1 

exp(-77TT) - Tl^ exp(-— ^). (32) 



V2^Ah 2A^' 5cV2^ 2A2j 
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The value of r.m.s. fluctuation at the horizon scale Ah can't reach value ~ 0.05 
which is necessary to achieve f3 ^ 10~® and correspondingly the cosmological 
mass fraction of primordial BHs ^ lO^'^ today. 

If, however, the spectrum has a peak at some scale, then primordial BHs 
are formed mostly in a narrow range of masses, near the mass that corresponds 
to this peak. The spectrum with a sharp maximum at some scale arises in 
some inflation models if the inflationary potential V{(j)) has a flat segment (see 
e. g. [1791 1180| ). As it was shown in |181j . the strong density fluctuations are 
generated in this models and collapse into SBHs. At the same time the spec- 
trum beyond the maximum may be of the standard Harrison-Zeldovich form and 
reproduce the usual scenario of large-scale structure formation in the galactic 
distribution. The primordial BHs production from the maximum may be con- 
nected to the formation of noncompact DM objects at the matter-dominated 
epoch pOO| . 

For the early formed primordial BHs, the process of mass increasing by 
the secondary accretion of DM is possible. Indeed, in the uniform Universe 
a primordial BH with mass inside the sphere containing total mass M 
produces the gravitational fluctuation ^oq = Mh/M at the epoch of matter- 
radiation equality t^q. Subsequent growth of this fluctuation obeys the well- 
known law for DM fluctuations. As a result the primordial BHs would be 
"dressed" by the DM halo with a steep density proflle, p oc r^^/^. We call 
this combined spherical volume of the "primordial BH -|-halo" by the "induced 
halo" . It was shown [48j that for primordial BH with a mass ~ IO^Mq the 
resultant mass of induced halo is of the order of mass of a typical dwarf galaxy, 
~ IO^Mq. During a hierarchical clustering of DM the induced halos enter into 
the haloes of ordinary galaxies. Many of these induced halos are capable to 
sink down to the galactic center during the Hubble time under the influence of 
dynamical friction. 

5.5 Clusters of BHs from closed domain walls 
5.5.1 General idea 

Some inflationary models suppose a creation of our Universe either near a max- 
imum of potential of inflaton field or near its saddle point (s) to realize a desired 
slow rolling providing a sufficient number of e-foldings (see details, e. g. in 
|201[ 1202] V As it will be shown below these models include the possibility of 
the formation of macroscopically large closed walls from a scalar field. After the 
end of inflation these closed walls collapse to BHs if these walls are large and 
heavy enough [1921 1193j . This mechanism is realized in well known models like 
the Hybrid Inflation [203j and the Natural Inflation 204 . A scalar field could 
be the inflaton itself or some additional fleld. 

First of all we consider a general mechanism of closed wall formation based 
on the quantum fluctuations near unstable point (s) like a saddle point or a 
maximum of potential of scalar field. An evolving scalar field may be split into 
a classical part, governed by the classical equation of motion, and the quantum 
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fluctuations [205j . To facilitate the analysis, let us approximate a potential near 
its maximum as 

V 



(33) 



m2 2 

where the maximum is assumed at = without the loss of generality. Then, 
a probability density to find a certain field value ((> has a form j206j (adapted to 
the considered case): 



(iP(0,r;0in,O) = d0 



7r(e2MT _ 1 



■ exp 



(34) 



Here a = pijcP' , = m? /iH and a = iJ^/2/27r, where the Hubble parameter 
~ V8^Fo/(3Mpi). 

Let us choose a positive value for the initial field, (f)^^ > 0, as illustrated in the 
Fig. [21 Then an average field value will increase with time, ultimately reaching 
a minimum of the potential at some value (p+ > 0. It means that a greater part 
of space will be finally filled with the field value = 0+. Meanwhile, the field in 
some (small) space domain could jump over the maximum due to the quantum 
fluctuations. In the following, an average value of the field representing this 
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fluctuation tends to an another minimum of the potential, < 0. As the 
result, the space at the final stage will be filled by vacuum (j>+ while some space 
domain is characterized by the field value (j) — (/)-< 0. If one starts to move 
from inside of the domain to the outside, the path would start from a space point 
with 0_ and finish at a space point with Hence, the path must contain a 
point with the maximum value of potential. It means that a wall is formed 
inevitably between such space domains and the "outer" space with (f> — 0+ 

[201 [Ml- 

The "dangerous" values of fluctuations are those with (p < 0. Such space 
domains will be surrounded by closed walls and if their number is too large 
it would strongly influence the dynamics of the early Universe. It is useful to 
calculate the probability of nucleation of these fluctuations. The latter depends 
on the initial field value 0in at the moment of creation of our Universe. The 
corresponding probability 

Po(0i„,r)= / dP{cj),T;cj)in,0) (35) 

to find the field value < at some space point for reasonable values of pa- 
rameters is represented in Fig. [3) This probability determines a ratio of spatial 
volumes with different signs of the field. This probability is highly sensitive to 
the initial field value ipi-a'. a closer to the potential maximum it is nucleated, a 
greater part of the Universe will be covered with walls at the final stage. 

If a fraction of space surrounded by the walls is not very large, the resulting 
massive BHs, which are formed from the walls, could explain the early formation 
of quasars [208j . In opposite case a contribution of BHs into the DM density 
is unacceptably large [207] . A further increasing of wall content leads to the 
wall-dominated Universe [209] . 

A mass and space distribution of resulting BHs strongly depends on a specific 
model and parameters of involved Lagrangian. It is instructive to demostrate 
the mechanism of massive primordial BHs production in the framework of the 
hybrid inflation model 203J following to the results of [207j . 

5.5.2 SBH formation in the Hybrid Inflation 

A potential of the hybrid inflation model according to [203) has the form 

V{x, a) = >c' (m' - ^) ' + ^X'a' + \m^a\ (36) 

Inflationary expansion of the Universe takes place during a slow rolling along 
the valley x = Oi <J > (Jc- When a fleld a passes through the critical point 

= V2-M, 
A 

a motion along the line x = 0, cr < CTc becomes unstable. As a result, a 
fleld X quickly moves to one of the minima, x± = ±2Af, cr = 0, see Fig. 2] 
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Figure 3: Part of space occupied by an another vacuum state, depending on 
the time measured by e-folds. The initial field value is (f)in = 0max + 3iJ, where 
^max ^ field value at a maximum of the potential. The parameter m is chosen 
to be m = 0.3-ff, where H is the Hubble parameter at the top of potential. 

where potential l|36p for the hybrid inflation model is represented. Inflation 
is finished with the intensive field fluctuations around one of an accidentally 
chosen minimum. 

It is clear now that this well elaborated picture suffers a serious problem 
nevertheless. During an inflationary stage, when the fields a and x move classi- 
cally along the line x = 0, the space is divided into many causally disconnected 
regions. The values of scalar fields within these regions are slightly different due 
to quantum fiuctuations. An each e-fold produces approximately ~ 20 space 
domains. Hence there are about e^*° ~ 10^* space domains right before the 
end of inflation. The values of field are chaotically distributed around the point 
X = 0, (7 = (7c- The domains with a field value x < tend to the left minimum 
X- = — 2A/, a — 0. Another part fall to the right minimum x+ = -I-2M, a — 0. 
A lot of walls between these domains appear and we come to a well known 
problem of the wall-dominated Universe |209j . 

The only way for our Universe to evolve into the modern state is to be 
created with an initial field value Xin ^ at the beginning of inflation. During 
inflation, the field x must slowly approach to the critical line x = 0. If, in the 
middle of inflation, an average field value approaches to the critical line x = 0) 
the field in some part of space domains could cross this line due to fluctuations. 
In the future, these domains will be filled by vacuum, say, x- surrounded by 
a sea of an another vacuum x+- These two vacua are separated by a closed 
wall as it was discussed above. A number of these walls strictly depends on the 
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Figure 4: A form of potential for the hybrid inflation model. Arrows indicate 
directions of classical motion of the fields. A black dot indicates the critical 
point. 
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initial conditions at the creation moment of our Universe. 

Let us estimate the energy and size of closed walls in the framework of the 
hybrid inflation. To this end, let us suppose that a field in some volume crossed 
the critical line at e-folds number N before the end of inflation. Its size is about 
the Hubble radius, , and it will be increased in times during inflation. A 
surface energy density of the domain wall after inflation defined by the potential 
dMl) is 

. = ^-^.M\ (37) 
Thus an energy i?waii contained in the wall after inflation is at least 

E^^n ^ 4vre {R-'e'^f = ^V2^e'^ . (38) 

A numerical value N varies in the interval (0 < iV < Nu ~ 60). These walls 
will collapse into the BHs with mass A/bh — -E'waii (for details see jl93] ). 

Let us estimate masses of these BHs for a representative values of parameters 
K = 10^2 and ^ 10^6 GcV. If iV = 40, we obtain for the mass of formed BH 

A/bh ^ 3 10^^ GeV - lOOM©. 

A similar estimation for a smallest BHs, which are created at the e-fold number 
N — 1 before the end of inflation, gives 

A/^BH, small — lO^A/pj. 

Therefore, the hybrid inflation leads to the formation of BHs with mass in a 
rather wide range, AfsH > 10^^ GeV with an upper limit of ^ IO^Mq. An 
amount of massive BHs depends on how close an average field value approaches 
to the critical line x = 0- The last, in turn, depends on the initial conditions 
and specific values of parameters of this model. 



5.5.3 Evolution of massive cluster of primordial BHs 

The main lesson of previous discussion is that the massive primordial BHs are 
formed after infiation as a rule rather than as an exception. As was shown in 
|194j the mechanism described above leads to the nucleation of groups (clusters) 
of primordial BHs rather than to a single primordial BH. Here we discuss the 
idea that these clusters of primordial BHs could be the seeds for early galaxy 
formation |208j . Our consideration is based on the results published in [1921 
119311206] . where the "Mexican hat" potential for the scalar field was considered. 
The calculations with the "Mexican hat" potential give the plausible parameters 
of clusters of primordial BHs. 

In numerical calculations, the cluster of primordial BHs is modelled by a 
spherically symmetric system with a radius r < ct, consisting of primordial BHs 
with a total mass Mh inside the radius r. To fix the cosmological evolution 
after infiation it is also supposed that the density of radiation is pr , the density 
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of ordinary DM is pdm and the density of A-term is pA- The main, most 
massive BH in a chosen cluster is in the center of the system. The density of 
radiation (and obviously the density of A-term) is homogenous. Therefore, the 
fluctuations induced by the primordial BHs are classified as entropy fluctuations. 
The scale under consideration is smaller than the horizon scale, and we use the 
Newtonian gravity with the prescription |210j of treating the gravitation of 
homogenous relativistic components, p — s- p + Zpjc^ . 

The evolution of any fictitious spherical shell of considered cluster with an 
initial radius ri < r obeys the equation 

<Pr _ G{Mh + A/dm) STrGp^r ^ ^irGpAr 

dt^^ r2 3 + 3 ' ^ ' 

with an approximate initial conditions at the moment t.^: r = iJr, r{ti) — ri. 
In obtaining equation (|39p. it was taken into account that + ipr = Se^., 
ea + 3pA = — 2eA- Numerical calculations on the basis of equation ([M]) with 
an initial conditions indicate that the radius of the shell increases initially, and 
gradually its expansion decelerates. At some instant, a speed of expansion 
diminishes to zero. The shell is separated from the cosmological expansion and 
starts to shrink. All components of nonrelativistic matter — the DM, primordial 
BHs and baryons follow the dynamics of the shell surrounding them. As a result, 
the solutions of (|39p for shells with different initial radii provide us the density 
distribution of the DM and primordial BHs in the evolving cluster. 

It is suitable to rewrite the main equation by using a new variable b{t). 

r{t) = ^a{t)h{t) (40) 

In this parametrization, ^ is a comoving length, a{t) is a scale factor of the 
Universe normalized to the present moment as a(io) = 1 and the function b{t) 
characterizes the deflection of the cosmological expansion from the Hubble law. 
The variable ^ is connected to the mass of DM inside the considered spherical 
volume (excluding total BHs mass) by the relation Mqm = (47r/3)pDM(io)?^) 
where pDM(io) is the modern DM density. The function a{t) obeys the Friedman 
equation, which can be rewritten as a/a = HqE{z), where z = a~^ — 1 is the 
redshift, Hq is the current value of the Hubble constant and function E{z) has 
a form 

E{z) - [flrflil + zf + an.0(l + zf + r!A,o]'/', (41) 

where ri,-^o is a density parameter of radiation, ^m,Q — 0.3, f^A.o — 0.7, and 
h = 0.7. Using the second Friedman equation (for a), one can rewrite p9p as 
follows 

'^^^ . '^^ ^ , f^+^h ,\ an.o(l + z) 



where function 



^^'^ = ' ^z 
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and a value of fluctuation 6h ~ Mh/Mj^M- Equation (|^ is equivalent to those 
obtained in [21 Ij in the case f^A — 0. We start tracing the evolution of cluster 
at a high redshift Zi when the considered shell crosses the horizon r ^ ct. 

Initial mass profile Mh(ri) of the cluster of primordial BHs was calculated 
following to formalism developed in [194] and has the form presented in the 
Fig. [5] (upper panel). This numerical result is a starting point for the following 
analysis. The mass Mi^Mifi) of DM inside the same spherical shells is also shown 
for a comparison. The radius is the physical size of sphere at the moment ti 
and the temperature Ti corresponds to the moment when the sphere is crossing 
the cosmological horizon. Note that the radiuses of shells in Fig. [5] are taken at 
different instants ti. Therefore, the shown mass of uniformly distributed DM not 
follows to the law Mdm c>c as it must be for a fixed moment, common for all 
spheres. Physical size at the temperature Ti is smaller than that in the modern 
epoch in To/Ti times, where To = 2.7 K. A density in the local center is so high 
that a lot of primordial BHs appear to be inside their common gravitational 
radius rg = 2GM/c^. According to numerical calculations with specific values 
of parameters, a total mass of primordial BHs appears inside the horizon from 
the beginning is 2.71O''M0. Hence, the massive BH with this mass is formed at 
first. 

A widespread prejudice is that a mass of formed primordial BH could not 
exceed a total mass under the horizon at the same time, M ^ {t/tp\)Mp\ |212| . 
As a result, the mass of primordial BH seems must be not larger than lO'^M©. 
It is true indeed if one considers the formation of primordial BHs density fluc- 
tuations. On the contrary, in our case the mass of the formed primordial BHs 
is determined by the area of the closed walls, formed and stretched during infla- 
tion. Size of wall could be greater its horizon size in orders of magnitude. BH 
is nucleated when the wall collapse under its horizon. 

The following dynamics of the cluster of primordial BHs together with DM 
component could finally lead to the formation of galaxies provided that the 
numerous collisions of galaxies are taken into account. Results of calculations 
made below are applicable both for an inner part of the cluster, composed mainly 
from primordial BHs and collapsing at the radiation dominated stage, and for 
an outer regions of the cluster, where DM is a main dynamical component. 
The outer regions are detached from the cosmological expansion at the matter 
dominated epoch. 

The most early epoch in our calculations corresponds to the formation of 
central BH with a mass 2.71O'*M0 described above. The temperature of the 
Universe at that time is ~ 16 MeV. We suppose that DM has been already 
decoupled from radiation at this temperature. For example, neutralino DM 
particles with mass 100 GeV and slepton with mass 1 TeV decouples at the 
temperature ~ 150 MeV |213| . well before the considered epoch. Therefore, 
the neutralino DM can be treated as moving freely under the influence of 
gravitational force only. As a result the clustering of two-component medium 
(BHs-l-DM) is described by the equation ((42|) from the very beginning. The 
same situation is realized for DM composed of the super-heavy particles with 
mass ~ 10^^ — 10^"* GeV, which probably never been in the kinetic cqui- 
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Figure 5: In the top panel the initial mass profile Mh{ri) of the cluster of pri- 
mordial BHs and the mass profile M-o^iTi) of DM are shown. In the bottom 
panel, the r.m.s. density perturbation at the moment teq of matter-radiation 
equality is shown. Two cases are compared: the total density fluctuations pro- 
duced in the presence of the cluster of primordial BHs and those produced in a 
standard way as remnants of the inflationary stage. 
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librium with radiation. In the opposite case (for some another DM particles 
candidates), the growth of fluctuations in DM medium is suppressed by the 
friction on radiation until the kinetic decoupling, while the BHs are cluster- 
ing. The super-heavy particles are preferable for our model in comparison with 
the neutralino because their annihilation cross-section is very small, oc rn~^. 
In this case the problems with a huge annihilation signal from the considered 
protogalaxies are absent. 

All scenarios imply a further merging of BHs during the multiple coalescences 
of host minihalos [2141 12151 12161 1217] accompanied also by an additional gas 
accretion. It is worth to estimate the probability to find a nowadays galaxy 
without the SBH in the framework of the scenario discussed above. Induced 
galaxies containing massive BHs and ordinary small protogalaxies without BHs 
have masses around Mdm = lO^M©, while nowadays galaxies are as massive 
as Mdm — lO^^M©. Every collision of an induced galaxy with an ordinary 
protogalaxy produces the next generation of protogalaxy with a massive BH 
in its center. So about 10^ collisions have happened up to now. Suppose that 
the amount of induced galaxies is about 0.1% in comparison with the ordinary 
ones. Upper limit of probability to find a nowadays galaxy without a SBH is 
0.9991° ^ 4.5 10"^ Hence, even a very small amount of induced galaxies is 
able to explain the observable abundance of e. g. the AGN. 

6 Conclusions and Discussions 

The existence of SBHs at the centers of galaxies looks like an unavoidable fact. 
The SBHs play an important role in the formation and global evolution of 
galaxies and the intergalactic medium. Nevertheless a valuable theory of SBHs 
formation is still absent. It appears not so easy to explain the whole set of 
observations which are biased and sometimes contradictive. The SBHs mass 
correlate with the mass of the host galaxies, though exceptions are known. For 
example, a suspected SBH was found recently in the dwarf galaxy VCC128 
|218| . The maximum accretion rate onto the SBHs, at redshift 2; ~ 0.7 almost 
coincides with a peak density of luminous infrared galaxies at z ~ 0.8. One 
could expect a quick growth of already existing SBHs. At the same time, the 
most massive BHs are found much earlier, at z > 6.0. 

At the recent epoch we are witnesses of impressive competition of two main 
directions in BH investigations. Several decades ago an idea was accepted that 
all BHs, if exist, result from the supernova explosions. Formation of galax- 
ies seemed definitely precede a formation of BHs. But gradually this picture 
becomes more and more tangled. New observational data during the last 20 
years indicate that almost all massive galaxies contain SBHs with a mass in the 
range 10^ — lO^M©. Moreover, the discovery of quasars at high redshifts means 
a combined evolution of the galactic nuclei and SBHs. The ubiquity of SBHs 
in galaxies may be the indication of some special and universal mechanism for 
SBH formation: simultaneously with the galaxies or even before the formation 
of galaxies. 
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The nowadays observations reveal a very early formation of SBHs: they 
already formed at the epoch of 800 million years after the Big Bang. Consid- 
erable efforts are needed to explain so early existence of SBHs originated from 
the remnants of stars. It seems that an idea of Population III stars with mass 
exceeding IO^Mq is able to solve the problem. Nevertheless, it is not excluded 
that astrophysical scenario for SBH formation will not succeed. 

Meantime new direction of primordial BHs research becomes more and more 
popular. The idea of primordial BHs is not very new — the first paper were 
pubhshed in 1966, [SHI [HZ]. It was supposed that primordial BHs were produced 
from density fluctuations long before the star formation. These BHs are very 
light, Mbh ^ 10^^ g, and should be quickly evaporated due to the Hawking 
radiation. These BHs could not merge effectively into the rather massive ones. 
Nevertheless, an idea of primordial BHs was survived. Several absolutely dif- 
ferent mechanisms of massive and even supermassive primordial BHs formation 
has been revealed and widely discussed [Ml Ea [HQl [Ha EHl HHH UMl fTM 

nmnmnm\Tmnm\mi\m[rm\m\ . 

Competition between astrophysical and cosmological schemes of SBHs origin 
is going on. New observational data are analyzed from the both points of view. 
Both astrophysical and cosmological scenarios posses as advantages and some 
defects. The astrophysical scenarios of SBHs formation has rather firm basis but 
suffer from the lack of time and shortage of matter supply for building the SBHs 
in galaxies by means of accretion and merging. Meanwhile, the cosmological 
scenarios are based only on the supposed mechanisms but provide the way of 
early SBHs formation. It is not excluded that both scenarios of SBHs formation 
were realized in the Universe. More elaborate and predictive models and new 
detailed observations are requested to resolve the SBH enigma. The upcoming 
gravitational wave laser interferometers provide a very promising new channel 
for the exploration of SBH problem. 
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